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Abstract 

The strange quark mass is determined from a study of the correlator 
(0|T ( J s Js) |0) (where J s = d^sj^u, the divergence of the strange vector 
current) using a family of finite energy sum rules recently shown to be very 
accurately satisfied in the isovector vector channel. It is shown that the match 
between the OPE and hadronic representations for the sum rules employed is 
extremely good once one fixes the overall scale, (m s — m u ) 2 , which normal- 
izes the OPE side. A value m s (4 GeV 2 ) = 115 ± 8 MeV, corresponding to 
m s (l GeV 2 ) = 159 ± 11 MeV, is obtained, with theoretical and experimental 
uncertainties contributing roughly equally to the combined error. 

12.15.Ff,11.55.Hx 
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I. INTRODUCTION 



The light quark masses, despite being fundamental parameters of the Standard Model, 
are surprisingly poorly known at present. Since, in the standard ("large condensate") sce- 
nario for the chiral counting in Chiral Perturbation Theory (ChPT), however, the ratio 
m s / (m u + rrid) = 24 .4 ± 1.5 is known fixing either m s or m u + is sufficient to deter- 
mine the other. In this paper we focus on m s , which has been the subject of several recent 
investigations, employing both QCD sum rules [P~|TI|| and lattice QCD fl2H23f. 

On the lattice, reliable results are presently restricted to the quenched approximation 
and vary somewhat depending on the fermion implementation and on whether one chooses 
itlk or rrifp when matching to physical meson masses. Unquenched calculations are at a more 
preliminary stage, but suggest unquenching is likely to lower quenched results 



Recent quenched results for m s (where here, as in what follows, all masses, unless otherwise 
specified, will be quoted at a scale fi = 2 GeV in the MS scheme) are 110 ± 20 ± 11 MeV 
Q, 122 ±20 MeV |3|, 95 ± 16 MeV @, 115 ±2 (143 ±6) MeV |7| (for m K (m ) input), 



111±12 MeV |18[, 98.1±2.4 MeV (update as quoted by Kenway @), 130±2±18 MeV 
[3g, 129 ±12 MeV [|22j and 95 ±26 MeV |23j. For dynamical (JV> = 2) simulations, the CP- 
PACS results m s ~ 70 (80) MeV (with m K (m^) as input) and SESAM result 151 ±30 MeV 
[16| do not appear compatible. (Note, however, that there are unresolved issues of chiral 
extrapolations, renormalization constants, and the lack of extrapolation to the continuum 
limit in the SESAM result; for a further discussion of these issues see Section 6 of the review 
by Bhattacharya and Gupta [12|, and Allton ||14|| .) 

Recent attempts to extract m s using QCD sum rules are of four types. Two employ flavor- 
breaking differences of vector and/or axial vector current correlators, either the isovector- 
hypercharge difference or the isovector-strange difference lT|| . The others involve 



analyses of the strange scalar |2],f§-§| and strange pseudoscalar g|| channels. 

For the vector current isovector-hypercharge difference, isospin-breaking corrections are 
crucial due to the high degree of cancellation on the hadronic side of the sum rule. These are, 
for natural reasons |7,24] , larger than one would naively expect, and can only be obtained 
indirectly 0. The high degree of cancellation also magnifies experimental and theoretical 
uncertainties, leading to significant errors on m s (m s = 107± ~ 40 MeV ||25|| , if one takes 
the isovector spectral function from hadronic r decay data P6],p7[). 



The isovector-strange difference, which employs r decay data for the isovector and strange 
spectral functions, is subject to complications resulting from (1) very poor convergence of 
the perturbative series for the longitudinal contributions on the OPE side |9|JT0|], and (2) 
strong cancellation, which again significantly magnifies the effect of the (otherwise rather 
small) experimental errors. Combining errors in quadrature and translating to the scale 2 
GeV, one finds m s = 163 ± 50 MeV [|, 143 ± 42 MeV and 157 ± 37 MeV |T 



The third approach involves the correlator (0|T (^Ja,s Ja,s) |0) (where Ja, s = d^s^^u) 
Hill. In this case, the K pole contribution to p(s) is known, but contributions from the 
continuum (i^(1460), i^(1830) resonance region) are not. Both analyses of this type employ 
a sum-of-Breit-Wigner modulation of the near-threshold tree-level ChPT form to deal with 
this problem, together with some assumptions about the relative size of the two unknown 
resonance decay constants. As discussed in Refs. P,p8| , this type of normalization procedure 
for a continuum spectral ansatz can involve potentially significant systematic errors. For 
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example, in the scalar strange channel, it leads to a rather poor representation of the actual 
spectral function (compare the model versions of Refs. H[| with that of Ref. H). In addition, 
because both analyses employ Borel transformed (SVZ) sum rules |[29|| , they involve the 
conventional "continuum" (or local duality) approximation to p(s) above s = Sq (where 
the "continuum threshold", So, is a parameter of the analysis). Since, for Borel masses, 
M, in the stability windows of both analyses, M 2 << is far from being satisfied, there 
are potential systematic errors resulting from violations of local duality in the "continuum" 
region. (That contributions from the "continuum" region are not numerically negligible can 
also be inferred, e.g., from the figures of Ref. ||). It is worth noting that violations of local 
duality are, in general, expected in kinematic regimes where typical resonance widths are 
comparable to or smaller than typical resonance separations. This has been quantified in 
the isovector vector channel, where the hadronic spectral function is known experimentally; 
violations of local duality in that channel, even at scales ~ m 2 , are seen to be, in general, 
rather large |3"0fl . Taking the results from Refs. which employ Aqcd — 380 MeV (which 
best corresponds to recent experimental values for a s (m 2 ) p6| , |27[| ), one finds, for m s , the 
values ~ 100 MeV (see Figure 5) and 115 ± 19 [Q (where the errors do not reflect the 
systematic uncertainties discussed above). 

The last approach involves the correlator, relevant to the strange scalar channel, 

IT(g 2 ) = i J dxe**(0\T(j a (x)Jl(0))\0) 

= i{m s -m u f f dxe iq - x (0\T(su(x)us(0)) |0) , (1) 



which has been analyzed, in Refs. using the Borel transformed sum rule method. 

This channel is a favorable one for sum rule analyses because the hadronic spectral function, 
p(s), turns out to be rather well-determinable (albeit indirectly) from experimental data. 
The reason is that since, experimentally, Kir scattering is known to be purely elastic up to 
s = (1.58 GeV) 2 pl|j3^| , the phase of the timelike scalar form factor, d(s), is known up 
to this point. d(s), moreover, satisfies an Omnes representation 



d(s) = d(0)p(s) exp 



s r°° , 5(t) 
1 dt- 



7T Jth t(t — S — it) 



(2) 



where p(s) is a polynomial normalized to 1 at s — 0, 5(s) is the phase of d(s), and the 
normalization d(0) is known from ChPT and K e ^ ||33|1 . In order to satisfy quark counting 
rules for the asymptotic behavior of d(s), 5 must — > (N + l)ir as s — > oo if p(s) is of degree 
N. In previous analyses p(s) has been assumed to be 1 (or, equivalently, very slowly varying 
in the range of s important to the analysis). The fact that, at the edge of the experimentally 
explored region, s = (1.7 GeV) 2 , the 1 = 1/2 Kir phase shift is very nearly n |3T| , |32| , |34| , 
the asymptotic value required if p(s) = 1, suggests, first, that this assumption is rather 
plausible and, second, that the treatment of the high energy behavior of the phase, S(s), 
employed in Ref. ||, i.e. 5(s) = ir for s > (1.7 GeV) 2 , is sensible. We will also see below 
that, making these assumptions, one obtains extremely well-satisfied sum rules, providing 
strong post facto justification as well. With these assumptions, the Kir contribution to p 
is completely determined, and multiparticle contributions begin only above s ~ 2.5 GeV 2 . 
The resulting Kir contribution to p(s) is, moreover, much smaller in the i^o(1950) than in 
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the Kq(1430) region (see Fig. 2 of Ref. 0). Since the Kq (1950) Ktc branching fraction is 
known to be ~ 50% [Q, the remaining multiparticle contributions to p, in the i^Q (1950) 
region, will also be small. This means that, to the extent that one may expect multiparticle 
contributions to p to be resonance dominated, (1) such contributions will play only a small 
role in the full spectral function out to s ~ 4 GeV 2 , and (2) in this region, the magnitude of 
these contributions can be obtained by combining information from the Omnes-generated 
Ktc component and the known experimental _K"q(1950) Ktc branching fraction. We do not 
quote results from the analyses of Refs. [Q,|J since the spectral ansatz employed there has 
been shown to be incompatible with the Omnes representation [J5J. Results from Refs. []5]|)j 
are m s = 91 — > 116 MeV || and 116 ± 22 MeV ||. The main potential uncertainty in 
these analyses is that associated with violations of local duality in the "continuum" region, 
s > sq. That continuum contributions are non- negligible may be seen from the lack of a 
stability plateau for m s in the analysis of Ref. ||. Finite energy sum rules constructed 
using the extracted values of m s from Refs. and the Omnes-generated spectral ansatz 
of Ref. || (also used in Ref. ||) as input have been shown to be moderately well satisfied 
||30|| . The residual deficiencies could arise from either shortcomings in the spectral ansatz or 
approximations on the OPE side of the sum rules. We will investigate this question in what 
follows. 



II. THE STRANGE SCALAR CHANNEL ANALYSIS 

Unitarity and analyticity ensure that correlators like n(g 2 ), defined above, satisfy ei- 
ther appropriately subtracted conventional dispersion relations or finite energy sum rules 
(FESR's). Borel transformation of the former leads to the conventional |29|] (SVZ) form of 
QCD sum rules. In what follows we choose to employ FESR's, which are based on Cauchy's 
theorem, rather than the Cauchy representation theorem, and work with the "Pac-man" 
contour, which traverses both sides of the physical cut, and is closed by a circle of radius so 
in the complex-s plane. One then has 

f s o —1 y 

/ ds p(s)w(s) = - — ; f dsw(s)U(s) , (3) 

J Sth ZTCl J\s\=sq 

valid for any function w(s) analytic in the region of the contour, where p(s) is the spectral 
function and s t h the relevant physical threshold. To make this relation useful, one wishes 
to use the physical spectral function, whose extraction was described in the Introduction, 
on the LHS, and the OPE for the correlator on the RHS. Since the OPE is proportional to 
(m s — m u ) 2 , this provides a method for extracting m s — m u . 

To make the results of this analysis reliable, Sq should satisfy two criteria: (1) it should 
be large enough to make the OPE representation of the RHS accurate, and (2) it should 
be small enough that one still has experimental information on p(s) for s < s on the LHS. 
The latter criterion, in the present case, restricts us to the region up to and including the 
_K"o(1950) (i.e., so < 4 GeV 2 or so). Since, even at s = 4 GeV 2 , resonance widths and 
separations are comparable, one expects local duality to be violated, creating potentially 
significant uncertainties on the OPE side. It has been argued, however, that even if local 
duality is not exactly satisfied for given timelike Sq, the OPE may, nonetheless, provide a 
good representation of the correlator on most of the circle \s\ = Sq (i.e., except possibly for 
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a region of hadronic scale near the timelike real axis) [p6[ . The isovector vector channel, 
for which the spectral function has been extracted from hadronic r decay data, provides 
explicit empirical confirmation of this argument. First, the conventional FESR treatment 
of these decays |37|] is very well satisfied, presumably because it involves a weight function 



(determined by kinematics) with a double zero at s = so which cuts out contributions from 
the potentially dangerous part of the contour. Second, one may demonstrate that using 
weights s k (with no such zero at s = so) leads to significant violation of local duality, even 
at scale Sq = m 2 , whereas FESR's based on weights having either a single or double zero at 
s = Sq are all extremely well satisfied, even at scales significantly below sq = m 2 [BUJ . Using 



such FESR's one may, thus, hope to construct sum rules in which the use of the "continuum" 
ansatz for the high energy portion of p(s) may be avoided. 

In what follows, therefore, we study FESR's for the correlator, II(g 2 ), using weights of 
the form 

■»« - i 1 - i) ( A - i) • (4) 

In Eq. (|j), A is a free parameter, to be chosen in such a way as to optimize the reliability 
of the analysis. Since (1) instanton effects are known exactly only in dilute gas single- 
instanton-background approximation [p8| , p9| and (2) the instanton density used to obtain 
these results is known to be altered by quark and gluon condensate effects |ItJ , the magnitude 



of the resulting instanton contributions (as stressed by the authors of Refs. |38| , |39| ) can be 
considered only a rough guide. For this reason, we restrict our attention to weights which, 
after integration, yield instanton effects (including corrections to the coefficients of both the 



unit operator and quark condensate in the OPE, as given in Ref. |3*9]|) which are small for all 



So employed. While any A > 1 turns out to satisfy this criterion, A = 2 creates an optimal 
suppression and we will display results below for this value of A. 

We now discuss the input on the hadronic side of the sum rules. The Ktt and multiparti- 
cle portions of p(s) are obtained, using the Omnes representation, experimental Ktt phases, 
and the experimental Ktt branching fractions of the Kq(143Q) and Kq(195Q), as described 
in the Introduction. We employ two slightly different implementations of this procedure, 
the difference in the resulting m s values serving to quantify the effect of uncertainties in 
our knowledge of the hadronic spectral function. In the first implementation, the fit to the 
experimental Ktt phases obtained in Ref. || is used as input to the Omnes representation, 
and multiparticle contributions are allowed in both the Kq(1430) and i£o(1950) regions, 
constrained by the Particle Data Group's Ktt branching fractions [R3]. In the second imple- 
mentation we employ the LASS fit |31"| , |32"| to the Ktt phase shifts and inelasticities. For this 
fit, the form of the parametrization of phase shift data used is slightly different from that of 
Ref. 0. In addition, since the fit to the Ktt amplitude is still purely elastic in the _K"q(1430) 
region, in employing this implementation we allow multiparticle contributions to the spec- 
tral function only from the _K"q(1950) region. (Note that the background effective range 



parameters, a, 6, and 1^(1430) resonance parameters, m, T, were mis-recorded in Ref. |JT 
and should actually be m = 1.412 GeV, T = .294 GeV ©gTJ, and a = 2.19 GeV~\ 
b = 3.74 GeV -1 [fffi| . It has been checked that these provide a satisfactory representation of 
the amplitude and phase shift data of Ref. |BI| 



On the OPE side, it is convenient to work with d 2 Il(Q 2 ) / (dQ 2 ) 2 = II" (Q 2 ), which satisfies 
a homogenous RG equation, so all scale dependence can be absorbed into the running mass 
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and coupling. To improve the convergence of the perturbative series, and significantly reduce 
residual scale dependence, we follow standard practice and evaluate the OPE side using 
the "contour improvement" prescription of Ref . , first summing logarithms through the 
renormalization scale choice fi 2 = Q 2 , and then performing the integral around the circle 
\s\ = sq numerically, using the exact solutions for the running mass and coupling which 



follow from the four-loop-truncated (3 and 7 functions [43,44 



The expressions for the dimension D = 0,2,4,6 contributions to U"(Q 2 ) are all known 
from previous work and, for /1 2 = Q 2 , are: (1) for D — 



n"(Q 2 



3 cm., 



D=0 



8tt 2 Q 2 
+a(Q 2 ) 



'i + H^!W) 2 



/5071 
V 144 



¥«3) 



, 3 / 1995097 



V 5184 



36 216 v ; 12 sw 



(5) 



where a(Q 2 ) = a s {Q 2 )/7T and C(-^0 is t ne Riemann zeta function; (2) for D = 2 
3(m s -m u ) 2 (Q 2 )m 2 (Q 2 ) 



n"(Q 2 ) 



D=2 



4tt 2 Q 4 



(6) 



(3) for D = 4 |,| 

n"(Q 2 



(m s -m M ) 2 (Q 2 ) f . . / 23 . ~,\ 1 / 121 . 
D=4 = 1 Q 7 { ~ ) {2{m s uu) (l + -a(Q 2 )) - -I G (l + —a(Q 2 , 



64 



7n [ 



a(q 2 



155 
~24~ 



(7) 



where Iq and I s are the RG-invariant D = 4 condensate combinations of Ref. and (4) 
for D = 6 H, 



n"(Q 2 ) 



(m s -m M ) 2 (Q 2 ) r _ 32 2 / 2 2 x 

' 3m s {gucrFu) — — 7T pvsA& \ \ ss ) + (wu) + 9{ss){uu)J 



Q* 



(8) 



where p^5A is a parameter describing the deviation of the four-quark condensates from their 
vacuum saturation values. In the above equations, factors of m u ^ have been dropped, except 
in the overall normalization (m s — m u ) 2 . 

For the input required on the OPE side of the sum rule we employ the following: 
a s (m 2 ) = .334 ± .022 |26j (needed as initial condition for the RG running of a s ); for 
the D = A condensates, (a s G 2 ) = 0.07 ± .01 GeV 4 g, (m u + m d ) (uu) = -/ 2 m 2 , and 
0.7 <r c = (ss)/(uu) < 1 for the D = 6 condensates, (gqaFq) = (0.8 ± 0.2 GeV 2 ) (qq) 
and Pvsa = 5±5 (i.e., allowing, to be conservative, up to an order of magnitude deviation 
from vacuum saturation). We evaluate m s — m u by optimizing the match of the OPE and 
hadronic sides of the sum rules in question in the range 3 GeV 2 < s < 4 GeV 2 , for which 
the dominant D = contour-improved series is well behaved. To convert the extracted 
values of m s — m u to those for m s , we take the central ChPT value for m s /m u [[[]]. 
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The errors on the extracted value of m s associated with those on the input OPE param- 
eters, together with their sources, are as follows (where we list only those cases where the 
error is greater than 0.1 MeV): due to a s (m%.), ±1.0 MeV; due to Pvsa, ±0.7 MeV; due to 
(a s G 2 ), ±0.2 MeV; due to r c , ±0.2 MeV. More significant, on the OPE side, is the error 
associated with the truncation of the D = contributions at four-loop order. If we assume 
continued geometric growth of the expansion coefficients, including this estimate for the next 
order term lowers m s by 2.8 MeV. In the final results below we have included this estimate 
of the five-loop order term in arriving at the quoted value, and assigned an error of twice 
the difference between the four-loop result and five-loop estimate. If added in quadrature, 
this source of error swamps all others on the OPE side. 

The sources and errors on m s associated with uncertainties on the hadronic side of the 
sum rules (where greater than 0.1 MeV) are as follows: due to the error on the A"q(1430) 
Ktt branching fraction ]35| (in the case of the first implementation described above, where 
multiparticle contributions are allowed in the A"q(1430) region), ±3.5 MeV; due to that on 
m K*(i43o), ±0.6 MeV; due to that on r^'*(i43o), ±3.8 MeV; and due to that on T K *( 1950 ), 
±0.4 MeV. There is also a 0.4 MeV difference between the extracted central values for 
the two different implementations of the hadronic spectral ansatz; this uncertainty is, how- 
ever, swamped by that associated with the treatment of multiparticle contributions in the 
-K"q(1430) region in the case of the first implementation. We treat the latter as a systematic 
error common to both hadronic implementations, even though it is, strictly speaking, not 
consistent to include it in the case of the second implementation, based on the LASS fit, 
which is purely elastic in the A"q(1430) region. 

Taking the average of the central values associated with the two different hadronic im- 
plementations, and adding errors in quadrature, we then arrive at our final result, 



where the first error is associated with the OPE side of the sum rule (and is completely 
dominated by our estimate of the error due to trunctation of the D = series) and the 
second with uncertainties on the hadronic side (where errors associated with the inelasticity 
in the ^(1430) region and with the uncertainty in the ^(1430) width are roughly equally 
important and, again, dominate all other sources). This result corresponds to 



In Figure 1 we display the match between the OPE and hadronic sides of the A = 2 
sum rule in the fit window, for the optimized value of m s . The agreement is obviously 
excellent, increasing our confidence in both the assumptions about the form of the Omnes 
representation of d(s), and the value of m s that results. 



We have determined m s from a FESR analysis in which the hadronic side of the sum rules 
employed is well determined by experimental data, modulo questions about the presence 
or absence of a polynomial factor in the Omnes representation of d(s) and the high energy 
behavior of the phase S(s). The excellent agreement between OPE and hadronic sides shown 




(9) 




(10) 



III. SUMMARY 
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in the Figure argues strongly, albeit post facto, for the reliability of the input assumptions 
p(s) = 1 and S(s) = tt for s > (1.7 GeV) 2 . We have attempted to make conservative 
estimates of the errors associated with uncertainties on the hadronic and OPE sides of the 
sum rules and conclude that m s (a GeV 2 ) = 115 ± 8 MeV, compatible with results from 
recent quenched lattice calculations. It is worth noting that systematic uncertainties in the 
Borel transformed sum rule approach associated with the use of the "continuum ansatz" for 
the high energy portion of the hadronic spectral function are potentially significant in size 
on the scale of the errors quoted above. In Ref. 0, for example, a variation of ±9 MeV in 
m s (l GeV 2 ) results from variations of ±0.5 GeV 2 in the continuum threshold parameter, 
so- Such uncertainties are, of course, absent in the FESR approach, and help to keep the 
final errors small. 
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FIG. 1. The OPE (Iope) and hadronic (I had) sides of the A = 2 FESR defined in the text. 
The solid line gives the OPE side, the dashed line the hadronic side. 
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